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THE  EFFECT  OF  ELECTRIC  FIELD  ON  THE  FREE  ENERGY  OF  SbSj^O*! 


I.  INTRODUCTION 

The  knowledge  of  the  Gibbs  free-energy  coefficients  is  necessary  when  one 
attempts  to  determine  the  thermodynamic  properties  of  a  material.  This  has 
been  demonstrated  by  Devonshire  [1,2]  and  later  by  Jona  and  Shirane  [3]  for 
ferroelectric  materials.  An  Interesting  feature  of  Jona  and  Shirane' s  deri¬ 
vation  is  that  the  free-energy  coefficients  may  be  determined  by  several 
methods.  One  can  measure  the  dielectric  constant  and  polarization  along  the 
c-axis  versus  temperature  in  a  range  about  the  Curie  point  or  monitor  just  the 
dielectric  constant  versus  temperature  as  a  function  of  direct-current  (dc) 
electric  fields.  Both  methods  yield  the  coefficients  of  the  ferroelectric 
terms  of  the  free  energy  and  suggest  an  interesting  way  to  test  the  consis¬ 
tency  of  the  experiments.  Another  salient  feature  of  the  electric-field 
method  is  that  it  can  also  be  used  to  determine  the  phase  transition  order 
(3). 


The  comparison  looks  even  more  enticing  when  applied  to  alloys.  The 
basic  Devonshire  approach  has  been  extended  by  the  author  to  SbS^_xOxI  alloys 
[4-6],  SbSI^_xBrx,  and  on  PbZrxTi^_x0-j  alloys  [7].  Several  authors  [8,9]  have 
investigated  the  influence  of  a  dc  electric  field  on  the  shift  of  the  Curie 
temperature  of  pure  SbSI,  but  no  similar  attempt  to  find  the  shift  of  the 
Curie  pressure  has  been  reported. 


The  purpose  of  this  article  is  to  present  an  approach  to  SbS^_xOxI  using 
an  electric  field.  The  free-energy  coefficients  of  SbS^_xOxI  alloys 
determined  by  using  both  methods  may  be  compared  in  order  to  verify  their 
equivalence  and  contribute  to  the  state  of  knowledge  on  SbS^_xOxI  alloys.  In 
Section  II  the  necessary  relations  are  developed  to  solve  for  the  free  energy 
of  an  alloy  via  field-dependent  dielectric  constant  measurements.  In  Section 
III  the  experiment  is  described;  i.e.,  the  measurements  of  dielectric  constant 
versus  temperature  and  pressure  under  the  applied  dc  electric  field.  The 
experimental  results  are  also  given  in  that  section.  The  comparison  and  the 
discussion  of  the  results  conclude  the  work  in  Section  IV. 

II.  THEORY 

The  Gibbs  free  energy  of  a  crystal  having  an  axial  polarization  (such  as 
BaTiOj ,  SbSI,  etc.)  is: 


^i  "  3  ^  ^i  3  ^  ^i  3  * 


(1) 


where  all  the  coefficients  are  temperature  dependent  and  ^  is  the 
polarization  along  the  c-axis.  The  index  i  denotes  the  parents  of  the  binary 
alloy  (i  ■  1,2)  or  the  alloy  itself  (i  =0).  The  alloy  is  assumed  to  have  the 
same  symmetry  as  its  parents  [4-7].  The  coefficients  aQ ,  ,  yq  will  be, 

therefore ,  concentration-dependent . 

Assuming  the  following  properly  weighted  relation  among  the  free  energy 
of  the  alloy  and  those  of  the  parents,  one  can  proceed  to  calculate  the 
relations  among  all  the  coefficients  involved: 


2 


Aq  =*  A^x  +  A2(l  -  x)  +  A^x(l  -  x)  .  (2) 

Here  x  is  the  concentration  of  the  hypothetical  oxygen  parent  (SbOI),  and 
A^  is  a  quadratic  interaction  term.  One  can  assume  next  [4-7]  the  following 
explicit  form  for  A3: 


2  1  2  4  4  2 

A3=  3P2  3  +  3P223  +  3F2  3  +  Pi  3P2  3^ 


in  which  a3 ,  g3,  and  Y3  are  the  coefficients  of  a  quadratic  Interpolation 
energy.  One  further  assumes  [4-7]  that  P^  =  =  P^,  which  yield  the  coveted 

relations  among  all  the  coefficients: 


«0  =  “j*  +  "  x)  +  a3x(l  -  x) 


00  =■  BjX  +  $2(1  -  x)  +  33x(1  -  x) 


Y0  =  YjX  +  ^2(1  -  x)  +  Y3x(1  -  x) 


As  discussed  in  References  4-7,  two  options  to  determine  the  coefficients  an. 


3q,  and  Yq  come  to  mind  immediately: 


a.  Deduce  Oq,  3q,  and  Yq  directly  from  measurement  to  which  the  alloy  is 
subjected  directly. 


b.  Calculate  the  coefficients  of  the  other  alloys  from  the  interpolation 
constants.  This,  however,  requites  the  complete  knowledge  of  the  parents 
and  one  alloy  in  order  to  calculate  the  characteristics  of  all  other 


alloys. 
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In  both  cases  the  necessary  measurements,  as  derived  from  Devonshire 
theory,  are  shown  in  Table  la  and  the  assumption  is  made  that  out  of  all  the 
discussed  coefficients,  only  the  a^'s  depend  on  temperature  [3]  and  are  in  the 
regular  Curie-Weiss  form: 


a 


i 


(5a) 


(5b) 
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where  and  Tq  ^  are  the  Curie  constant  and  the  Curie-Weiss  temperature  of 
the  constituent  i. 


Jona  and  Shirane  [3]  considered  also  a  different  approach  that  leans 
heavily  on  the  effects  of  electric  field  on  the  dielectric  constant.  While 
is  determined  from  Eq.  (5)  in  the  familiar  routine  (see  Table  la),  and 
are  determined  from  the  following  relations: 


Ei,33(T)  "  Ei,33(T) 


12lei°33(T)]\ 

3  2 

(*»)V 


(6) 


dT 


c  ,i 


2Y. 


v2 


dE 


4?  <- 


(7) 


Here  eiE33(T)  is  the  dielectric  constant  at  a  nonvanishing  dc  electric  field  E 
while  ei^33(T)  is  its  value  for  E  ■  0,  both  for  the  ic^  constituent  and  the 
c-axis.  The  corresponding  gradient  of  the  Curie  temperature  with  electric 
field  is  dTCjl/dE.  Equations  (5)  through  (7)  are  sufficient  to  determine  o^, 
fjj,  and  Y^  and  are  summarized  in  Table  Ih.  A  bonus  originating  from  Eq.  (7) 
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Table  la.  The  free-energy  coefficients  and  the  corresponding  measurements  required 
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**Values  for  this  coefficient  are  given  at  some  arbitrary  temperature  Tr  for  the 
order  case.  Normally  Tr  is  chosen  as  the  room  temperature  (25°C)  but  in  this 
20°C  below  Trt. 


Is  the  potential  to  determine  the  order  of  the  transition.  While  Y^  is 
considered  always  positive,  3^  is  negative  for  a  first-order  transition  and 
positive  for  a  second-order  one.  A  change  in  the  sign  of  dT  H/dE  means  a 
second-order  transition.  As  this  quantity  approaches  zero,  it  denotes  a 
tricritical  point.  One  should  bear  in  mind  that  Eq.  (7)  is  correct  only  for 
first-order  transitions,  and  another  expression  should  be  worked  out  for  the 
second  order  [3] . 

When  Eqs.  (6)  and  (7)  are  written  for  an  alloy  using  Eqs.  (4)  to 
determine  the  interpolation  constants,  one  gets: 

E  0  12l£i°33(T)l4 

ei  ,33<T)  "  £i,33(T)  - - 3  -  t8ix  +  82(1  “  X)  +  83x(1  '  x)] 


and 

dT  .  21/2C.  Y,x  +  Y,(l  -  x)  +  Y„x(l  -  x)  1/2 

c  »*■  a  _ i.  r _  _* _ 2 _ 2 _ i  (<j) 

dE  4*  1  3  x  +  32(1  -  x)  +  S3x(l  -  x)J  '  w 

As  one  can  clearly  envisage,  Eqs.  (6)  and  (7)  are  typical  of  attaining  Oq,  3 
and  Yq  from  measurements  while  Eqs.  (8)  and  (9)  allow  one  to  calculate  the 
coefficients  of  other  alloys. 

The  extension  of  the  above-mentioned  considerations  to  include  nonvan¬ 
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ishing  pressures  is  straightforward  when  one  incorporates  the  Samara  [10] 
transformation  into  the  development.  The  pressure-dependent  coefficients 
starred.  It  follows  from  the  former  assumptions  [4]  that: 


(13) 


=  <T  -  T  )/ (dT  /dp) 
0,1  o,i 


«  ^/(dT 


o,i 


/dp)  . 


(14) 


Here  and  pQ  ^  are  the  Curie  constant  for  nonvanishing  pressures  and  the 
Curle-Welss  pressure.  As  implied  by  Eqs.  (13)  and  (14),  they  are  derived  from 
their  zero-pressure  counterparts  through  division  by  the  Ehrenfest  temperature 
gradient. 

The  derivation  of  expressions  analogous  to  Eqs.  (4),  (8),  and  (9)  is 
straightforward  and  will  not  be  included. 


III.  EXPERIMENT  AND  RESULTS. 

A  two-fold  experiment  was  conducted  that: 

a.  determined  the  polarization  and  dielectric  constant  as  a  function  of 


temperature  and  pressure. 


b.  measured  the  dielectric  constant  in  the  presence  of  an  electric  dc 
bias  field  as  a  function  of  temperature  and  pressure. 


The  first  part  of  this  experiment  has  already  been  described  elsewhere  [4,5] 
and  will  not  be  presented  in  detail  here. 

The  samples  used  in  both  experiments  were  grown  by  the  Material  Research 
Laboratory  of  the  Pennsylvania  State  University.  They  were  received  as  9-mm 
to  12-mm  diams  and  up  to  25-mm  long.  Samples  were  cut  using  a  diamond  saw  and 
polished  to  //600  before  electroding.  To  increase  their  mechanical  strength, 
they  were  impregnated  in  SPURR's  low-viscosity  epoxy  under  vacuum.  Then  the 
final  samples  were  reshaped  to  their  original  diameter  with  the  typical  final 
thickness  between  2  mm  and  6  mm.  Dynalloy  350  fast-curing  silver  paint  was 
used  for  the  electrodes.  All  samples  were  polarized  at  1.5-kV/cm  at  -11°C. 


The  dielectric  constant  was  measured  by  a  General  Radio  1683  automatic 
RLC  bridge  at  1  kHz.  The  samples  were  immersed  in  castor  oil.  Temperature 
was  controlled  between  -12  and  +60°C,  and  pressures  to  100  MPa  were  obtained 
using  a  pressure  pump  and  a  BLH  1200  digital  strain  indicator.  A  dc  field  was 
applied  by  a  filtered  power  supply. 

The  concentrations  of  the  oxygen  parent  used  (in  molar  percents)  were  0, 
4,  8,  12,  and  20.  As  the  samples  were  highly  conducting  neat  their  Curie 
temperature,  even  fields  of  10  V/cm  could  cause  breakdown.  Although  such 
fields  are  by  no  means  extreme,  some  breakdowns  did  occur.  This  prevented 
comparison  with  the  results  of  option  a  in  Section  II  of  the  experiment  for 
x  »  4,  8,  and  20%.  The  data  finally  used  were  measurements  on  pure  SbSI , 

SbSg .92^0 .08* »  an<*  SbSQ.880o.l21  t0  coraPare  the  Curie  temperature  and  pressure 
shifts  and  the  free-energy  coefficients.  The  data  were  plotted  as  the  Curie 
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temperature  and  pressure  versus  field  in  Figs.  1  and  2.  Figure  3  is  an 
example  of  the  actual  1/e^  33  measurements  as  a  function  of  E.  Linear 
regressions  of  1/e^  33  versus  T  (or  p)  were  calculated  (see  Fig.  3)  both  for 
the  last  three  points  below  Tc  ^  and  for  all  of  those  above  Tc  The 
straight  lines  obtained  were  then  intersected  to  find  the  Curie  temperature 
(or  pressure)  at  the  relevant  field. 


*  SbS0.92°0.08I 

■  PURE  SbSI 

■  SbS0  88°O.I2I 


-  Curie  temperatures  for  various  sbSL_xOxI  alloys  as  functions 
of  a  dielectric  constant  electric  fiela 


-  Curie  pressures  for  various  SbSt_xOxI  alloys  as  functions 
of  a  dielectric  constant  electric  field 


reel 


3  -  Reciprocal  dielectric  constant  vs  temperature  at  various 
dielectric  constant  fields  for  SbSg .8800.12* 
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IV.  DISCUSSION 

It  is  interesting  to  note  that  Tc  ^  changes  linearly  with  E  (see  Figs.  1 
and  2)  as  predicted  by  Eq.  (9)  for  each  concentration.  The  values  for 
dTc^/dE  are  also  in  fair  agreement  with  those  of  Toyoda  et  al.  [9].  It  is 
worthwhile  to  examine  the  dependence  of  the  Curie  temperature  gradient  with 
pressure  on  the  electric  field.  One  may  use  the  mathematical  identity: 


3T 


3E 


3E 


“  3p, 


)  ( 


3p 


c,i 


3T 


=  -1 


c,i 


(15) 


to  examine  this  dependence.  This  test,  when  applied  to  small  dc  electric 
fields,  can  check  the  consistency  of  the  experimental  method:  The 

l (dTc  ,i)/(3pc>i)]E  values  are: 

7.9  ±  2.0°C/GPa  for  pure  SbSI 
-4.3  A  1 .0°C/GPa  for  SbSg>q20(}.08^' 

-7.7  ±  2.0°C/GPa  for  SbS0,88°0.12I* 

The  value  for  pure  SbSI  is  in  fair  agreement  with  the  result  of  -4.7  A  2°C/GPa 
quoted  by  Samara  [11]  for  epoxy-bonded  samples. 

The  gradients  of  the  Curie  temperatures  and  pressures  with  electric  field 
for  each  alloy  are  summarized  in  Fig.  4  in  two  frames.  The  gradients  of  the 
Curie  tempertures  and  pressures  are  plotted  as  a  function  of  oxygen 
concentration  in  the  alloy.  The  data  presented  in  Section  III  are  plotted  as 
crosses  and  should  be  compared  to  the  values  calculated  from  Eq.  (7)  which  are 
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0.05 


0.03 

dpc,i  /MPo  -m 


Fig.  4  -  Gradients  of  the  Curie  temperatures  (A)  and  pressures  (B) 
as  functions  of  the  oxygen  concentration  In  SbSi_xO  I. 
Crosses  are  direct  measurements  (see  Section  III),  full 
circles  are  values  calculated  from  Eq.  (7)  using  the  free- 
energy  coefficients  at  zero  field  for  each  sample.  The 
curve  represents  the  result  of  the  Interpolation  scheme 
[see  Eq.  (9)]  when  applied  to  the  range  0<x<D.2. 


displayed  as  full  circles.  The  solid  curve  joining  the  triangles  represents 
the  dependency  predicted  by  the  interpolation  scheme  of  Eq.  (9)  when  applied 
to  the  range  0  x  <^0.2.  Due  to  the  shortage  of  experimental  data  with 
field,  it  is  difficult  to  test  the  validity  of  the  interpolation  method  when 
applied  to  the  field-oriented  results.  One  can  carefully  state,  though,  that 
the  interpolation  method  will  deliver  results  of  the  same  quality  regardless 
of  the  presence  of  field.  This  is  corroborated  by  the  proximity  of  the 
experimental  points  obtained  by  both  methods  (see  again  Fig.  4). 
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It  looks  in  Fig.  (4)  as  if  the  minimum  in  the  values  of  both  dTc  j/dE  and 
dp.  . /dE  may  indicate  that  the  alloy  system  possesses  a  tricritical  point 
phase  transition.  However,  the  evidence  is  only  circumstantial  at  this 
time.  It  is  noted  that  in  Fig.  4a  a  minimum  appeared  at  x  ■  0.08,  whereas  for 
x  -  0.04  the  alloy  SbSg  ^Oq  q^I  exhibited  the  largest  Tc  ^  and  pc  ^  (4]. 

This  situatin  may  be  coincidental.  The  magnitude  of  Tc  ^  and  pc  ^  depends  on 
the  coefficient  of  the  free  energy  (Eqs.  5  and  12),  and  dTc  ^/dE  and 
dpc  j/dE  upon  8^  and  coefficients  [Eqs.  (1)  and  (7)].  Therefore,  one 
cannot  correlate  the  peak  in  Tc  ^  and  pc  ^  with  the  valley  in  dTc  ^/dE  and 
dp.  j/dE,  respectively. 

The  high  conductivity  of  SbS^_xOxI  close  to  the  transition  point  proved 
to  be  detrimental  to  some  of  the  samples  used.  The  conductivity  depends,  to  a 
certain  extent,  on  the  growth  technique  and  the  choice  of  the  raw  materials. 
Some  of  the  points  raised  in  this  discussion  could  have  been  better  examined 
had  the  samples  been  less  conductive. 

One  should  consider  now  the  material  problems  contributing  to  the 
observed  errors.  The  orientation  of  the  samples  is  not  perfect  since  they 
were  grown  by  using  the  Bridgman  method.  Some  spread  in  results  is  expected 
because  of  the  misalignment,  and  one  can  expect  this  spread  to  broaden  if  the 
growth  mechanism  is  affected  by  the  addition  of  variou-  oxygen  quantities. 
While  the  pure  SbSI  and  the  SbSg^gOQ^i  samples  are  highly  aligned, 
SI)S0.9200.08I  an<*  s^s0.88°0. 12l  are  less  wel1  oriented.  A  sample  of 
SkS>0.96®0.04*  has  also  been  tested,  but  the  results  were  not  presented  due  to 
the  mediocre  orientation.  The  criterion  for  presenting  the  results  was  that 
the  worst  misorienttion  ot  needles  would  be  less  than  35°. 


i 
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One  should  also  consider  the  Instrumentation  precision.  The  accuracy  of 
the  capacitance  bridge  Is  better  than  0.1%,  and  the  field  calibration  was 
accurate  to  within  2%.  In  spite  of  these,  it  was  felt  that  the  accuracy  of 
the  transition  point  is  less  (mainly  due  to  the  ±0.2°C  inaccuracy  in  the 
temperature  and  ±0.01  MPa  in  the  pressure  readings).  The  transition,  although 
of  first  order  [12],  is  smeared.  The  considerable  slope  of  the  reciprocal 
dielectric  constant  with  temperature  (or  pressure)  contributes  to  the 
accumulated  error.  The  accuracy  of  the  measurement  was,  therefore,  secondary 
in  nature  to  this  effect  and  was  considered  adequate  for  the  purpose  of  this 
work. 
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